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I Abstract. We review the main results from line-driven (LD) models of winds 

<*' , from accretion disks in cataclysmic variables (CVs). We consider LD disk wind 

O ■ models in the hydrodynamic (HD) and magnetohydrodynamic (MHD limits. 

We discuss the basic physical conditions needed for a disk wind to exist and 
the conditions for the wind to be steady or unsteady. We also discuss how 
the line-driven (LD) wind structures revealed in numerical simulations relate to 
observations. In particular, we present synthetic line profiles predicted by the 
^ . LD wind models and compare them with observations. Our main conclusion is 

' that, despite some problems, line-driving alone is the most plausible mechanism 

I for driving the CV winds. This conclusion is related to two facts: 1) LD wind 

PsJ . models are likely the best studied wind models, analytically and numerically, 

T-H ' and 2) it is most likely that the predictive power of the LD wind models is much 

^-H ^ higher than of any other wind model so far. Preliminary results from LD-MHD 

■ wind models confirm that magnetic driving is likely an important element of 
' the wind dynamics. However, magnetic driving does not seem to be necessary 

^ I to produce a CV wind. The most important issues which need to be addressed 

by future dynamical models, regardless of driving mechanism, are the effects 
I ' of the position-dependent photoionization and the dynamical effects in three 

? , dimensions. 

a . 

^> : 

^ I 1. Introduction 

H ; 

I Winds in cataclysmic variables (CVs) exemplify very well a phenomenon of ac- 

cretion disks being accompanied by mass outflows. Other systems where this 
phenomenon occurs are active galactic nuclei (AGN) and young stellar objects 
(YSOs). In the case of CVs, key evidence for outflows comes from P-Cygni 
profiles of strong UV lines such as C IVA1549. However, the evidenc e for the 
winds is not limited just to the strong UV lines iLong Sz Kniggd (2003)- Under- 



standing the winds in CVs is important on its own right and because they have 
been the best observed outflows from compact objects and promise to provide 
us with insights into all disk outflows. The interpretation of data is usually 
based on fitti ng observed profiles to synthetic profiles calculat ed from kinematic 
models (e.g.. iMauche fc Ravmondlll987l: lDrewlll987l : fShlosma n fc Vitelld 1199^ : 
iKnigge. Woods, fc Drewlll99,llLong fc Kniggell2nn2r ) . We refer a reader to Fron- 
ing (2004, this volume) for a review of the observations and interpretations of 
CV winds. 

Magnetic fields, the radiation force and thermal expansion have been sug- 
gested as mechanisms that can drive disk winds. These three mechanisms have 
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been studied extensively using analytic as well as numerical methods. As a re- 
sult of these studies, theoretical models have been developed that allow us to 
estimate under what physical conditions each of these mechanisms is efficient 
in launching, accelerating and collimating disk outflows. Here we will focus on 
models of radiation-driven disk winds (section 2) and a hybrid model in which 
both radiation and magnetic driving is considered (section 3). 



LD HD Models 



More than three decades of studies of winds in hot stars provide us with a very 

good understanding of how line-dr iving produces powerful high velocitv winds 

('e.g.. lCastor. Abbott fc Kleinll97,l hereafter CAK: lFriend fc Abbottll"98^ : lPauldrach. Puis, fc Kudritzkill 



The key element of the ICAKi model is that the momentum is extracted most 
efhciently from the radiation field via line opacity. The Eddington limit, L^dd = 
4TTcGM/ae, is the maximum luminosity a spherical object of mass M may 
achieve before the radiation pressure mediated by photons scattering off free 
electrons becomes so large as to drive off the object's atmosphere and envelope. 
It is commonly the case that the effective cross section for photon scattering is 
greatly increased by the presence of Doppler-shifting bound-bound transitions. 
ICAK showed that the radiation force due to lines, F^"''^'^ can be stronger than 
the radiation force due to electron-scattering, F™'^''^ by up to several orders of 
magnitude (i.e., i?^ad,«^^rad,e ^ ^^^^ _ 2000). Thus even a star that radiates 
at around 0.05% (i.e., l/Mmax) of its Eddington limit can have a strong wind. 

Early models of radiation-driven disk winds have applied assumptions that 
either restrict the flow geometry or require the flow to be time-independent [e.g. 
IVitello k ShTosman (J.988) and a model for AGN wind by >Iurrav et al. (1995)]. 
Recently numerical models of 2.5-D, time-dependent radiation driven disk winds 
have been constructed for app lication to , in the first instance, CV d isk winds 



fe.g.. iPerevra. Kallman. &: Blondin 199^ : IProga. Stone. &: Drewl|l99^ . hereafter 
PSD 9 8: iProga Stone, fc Dr^ll999l . hereafter PSD 99: lProgJl999l ). In par 



ticular, IPSD 98l adopted numerical techniques to integrate the coupled HD and 
radiation transfer equations in order to study the multidimensional and time- 
dependent character of line-driven (LD) disk winds fro m first principles. For the 
HD equations, they used the well-tested ZEUS code ( Stone NormanI ll992D 



extended by the addition of a term in the equation of motion which accounts 
for the radiation force due to spectral lines of the form 

prad,l ^ f f ^edr\ ^^^^^ 



/E V C 

The term in brackets is the electron scattering radiation force (^F™'^'^) and M, 
the force multiplier, is the increase in opacity due to spectral lines. The integra- 
tion is over all visible radiating surfaces (S). Note that dJ- contains the total 
frequency-integrated intensity emitted at a given location. PSD 98 adopted the 
simple form for M which still underpins much modeling of OB star winds, i.e. 
M = kt~°', where t is proportional to the lo cal de nsitv divided by the local 
velocity gradient, and k and a are constants The maximum values of 

M, Mjnax determines the actual luminosity limit (the effective Eddington limit, 
LE/Minax) for which the radiation pressure mediated by photons becomes large 
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enough to drive off the object's atmosphere and envelope. Note that integration 
over angle of a nonlinear function of the velocity gradient tensor and the radia- 
tion flux is required to evaluate equation (1): Proea. Stone. Sz Drewl took great 
care in the numerical evaluat ion of these integrals by adopting angle-adaptive 
quadrature methods. IPSD 98l 's formalism allows the radiation from the central 
accreting star to be included both as a direct contributor to the radiation force 
and as an indirect component via disk irradiation and reemission. To spatially 
resolve the flow, we used a non-uniform (up to 200 x 200) grid in which the 
subsonic acceleration zone near the disk or stellar surface is well sampled. 

The primary outcome of the numerical studies is the confirmation that line 
driv ing can p roduce a supersonic, biconical wind from an accretion disk in CVs. 

IPSD 98l explored the impact upon the mass-loss rate, Mw and outflow ge- 
ometry caused by varying the system luminosity and the radiation field geom- 
etry. In their study, the system luminosity, L was defined as the sum of the 
disk luminosity and the central star luminosity, Ld and L^,, respectively (i.e., 
L = Ld + = (1 + x)Ld; where x = L^^/Ld). A striking outcome was that 
winds driven from, and illuminated solely by, an accretion disk yield complex, 
unsteady outflow (see the top panels in Figure 1). In this case, time-independent 
quantities can be determined only after averaging over several flow timescales. 
On the other hand, if winds are illuminated by radiation mainly from the central 
object, then the disk yields steady outflow (see the bottom panels on Figure 1). 
IPSD 98l also found that Mw is a strong function of the total luminosity, while 
the outflow geometry is determined by the geometry of the radiation field. For 
high system luminosities, the disk mass-loss rate scales with the luminosity in 
a way similar to stellar mass loss (e.g., compare the crosses and the solid-line 
curve on Figure 2). As the system luminosity decreases below a critical value 
(the Eddington factor, F = L/LEdd about twice 1/Mmax) the mass-loss rate 
decreases quickly to zero. 

The simulations also showed that regardless of the radiation geometry, the 
two-dimensional structure of the wind consists of a dense, slow outflow that is 
bounded on the polar side by a high-velocity stream (respectively 'slow wind' 
and 'fast stream', for short; see e.g.. Figure 1). Matter is fed into the fast stream 
from within a few central object radii. In other words, the mass-loss rate per 
unit area decreases sharply with radius. The terminal velocity of the stream is 
similar to that of the terminal velocity of a corresponding spherical stellar wind, 
i.e., Voo ~ a few Vesc, where Vesc, is the escape velocity from the photosphere. 
Thus the difference in geometry changes the wind geometry and time behavior 
but has less effect on and Voo- 

IPSD 9fll applied a more accurate treatment of the line force than IPSD 981 
by including all the terms in the velocity gradient tensor. Qualitative features 
of the new models are very similar to those calculated byj^SD 98. In particular, 
the more a ccurate calculations showed that models which displayed unsteady 
behavior in IPSD 981 are also unsteady with the new method, and gross properties 
of the winds, such as mass- loss rate and characteristic velocity are not changed 
by the more accurate approach. The largest change caused by the new method 
is in the disk-wind opening angle: winds driven only by the disk radiation are 
more polar with the new method while winds driven by the disk and central 
object radiation are typically more equatorial. 
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Figure 1. Maps o f poloidal velocity for a range of LD disk wind models 
(Figure 2 in IPSD 991) . The top panels are both models with a; = but with 
Ma = 10"^ M© yr~i (the left hand side panel) and Ma = tt x 10~^ M© yr"^ 
(the right hand side panel). The bottom two panels are results for models 
both with Ma = tt x 10^* Mq yr~^, but with x = \ (the left hand side panel) 
and X = 3 (the right hand side panel). The top two panels show the effect 
on the outflow geometry of increasing the disk luminosity alone, while the 
top right and bottom two panels show the effect of adding in an increasingly 
larger stellar component (a; = 0, 1 and 3) to the radiation field. Adding in 
an increasingly large stellar component causes the outflow to become more 
equatorial. Note that we suppress velocity vectors in regions of very low 
density (i.e., p less than 10~^° g cm~^). 



2.1. Nature of Unsteady Outflow 

The fact that the unsteady behavior observed in lPSD 9^ models has not changed 
with a more accurate treatment of the radiation force indicates it is indeed a 
robust property of line driven winds from disks. Why does increasing the radial 
com ponent o f the radiation force 'organize' the wind into a steady state? 

iPSD presented the following explanation; let r and z define position 
along a streamline in the wind in cylindrical coordinates. An increase of the 
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vertical component of the gravity, gz oc —zjir^ + z^)^/^ with height at a fixed 
radius r is the main driver of the unsteady flow. (We add that for a disk wind 
driven by the disk radiation is vertical and its density can decrease downstream 
only due to an increase in the velocity as there is no geometrical dilution.) The 
increase of the gravity with height can be significantly reduced if the streamlines 
are directed outwards from purely vertical (i.e., r increases with z). At the same 
time, this tilt also brings into play an increase of the horizontal effective gravity, 
g'r, along each streamline: gr oc [rf/r^— r/(r^ + z^)^/^], where rf is the radius on a 
Keplerian disk at which a streamline originates. However the increase of gr with 
r' is slower than the increase of g^, with z' because of the decaying centrifugal 
term. In other words, the line force can more easily maintain domination over 
gravity if the flow climbs the gentler gravitational hill in the horizontal direction 
as compared with the vertical direction. Furthermore, driving material along 
streamlines outward from the vertical causes density to decline with radius due 
to geometrical dilution alone as with p ~ 0.5 — 1.5 — this, very usefully, 
tends towards increasing the line force, thereby facilitating a better match with 

trends in gravity. 

We note that iProea. Stone. &: Drewl 's discovery of unsteady flow in models 
for X = has not been confirmed by Pereyra and his collaborators in their numer- 
ical s imulations ( Perevra. Kallman. &: Blondin lll997l . 120001 : iPerevra fc: Kallman I 
l20n.lD . Recentlv- IFerevra et al. I (l2004l 'l presented mathematically simple mod- 
els and used them to argue against the above explanation of unsteady na- 
ture of LD disk wind for x = 0. In particular, Pereyra et al. claimed that 
a gravitational force initially increasing along the wind streamline, which is 
character istic of disk wi nds, does not imply an unsteady wind. However, we 
find that iPerevra et al. I omitted the fact that the line force and consequently 
wind streamlines are coupled to the gravitational force and radiation field. 
We also note that Ivitello & Shlosman (1988) studied analytically LD driven 
disk wind model where streamli nes were vertical nea r the disk and diverging 
to spherical at large distances. IVitello &: Shlosmanl found that for a steady 
state solution to exis t it is n ecess ary to en force a radiation force term to in- 
crease with height IPSD 98l and IPSD 99l 's results are then consistent with 
IVitello &: Shlosmanl 's result because for the case where the wind is vertical 
near the disk only unsteady solution exist ((PSD 9& computed the radiation 
force self-con sistently and did not enforce any spec ial increase of it with height 
[See however IPerevra. Kallman. Blondin I ( 19971 )]. In summary, we conc lude 
Perevra et al. I 's mathematically simple models, if anything support lPSD 98l and 



PSD 99l 's interpretation of unsteady winds: when one considers the radiation and 



flow geometry appropriate for the x = case (i.e., radiatio n and flow are v ertical 
near the disk) then analytic analysis similar to that used bv lPerevra et al."1 shows 
that there is no critical poi nt fusing the ter r ninolo gy from analytic models of LD 
stellar winds as adopted bv lPerevra et al. ](|200l)). However, significant central 
radiation (as in the x > cases) makes a dynamical and geometrical change in 
the wind solution i.e., increases the wind inclination angle and a critical point 
can exist. 
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2.2. Testing Models Against Observations 

Applying IPSD 98l and IPSD 99l 's dynamical models to CVs one finds that ra- 
diation driving can produce disk winds consistent with the following observed 
properties of CV winds: (i) the flow is biconical rather than equatorial as re- 
quired by the absence of blueshifted line absorption from the spectra of eclipsing 
high-state CV, (ii) the wind terminal velocity is comparable to the escape veloc- 
ity from the surface of the white dwarf (WD), and (iii) the spectral signatures 
of mass loss show a sharp cut-off as the total luminosity in dwarf novae declines 
away from maximum light through the regime theoretically identified as likely 
to be critical. 

Additionally, LD disk wind models may explain the highly unsteady and con- 
tinuously variab le nature of the sup ersonic outflows in the NL binaries BZ Cam 
and V603 Aql (|Prinia et alJ l2nnn a''?l . The presence of a slow, dense transi- 
tion region between disk photosphe re and outflow in V347 Pup and UX UMa 
( Shlosman. Vitello. fc Mauchelll99(il : Ik nigge fc Drewl[l997l 'l may also be accom- 



modated within these same models. 

Building upon these numerical models, one can compare the results of these 
mode ls with the a nalvtic results readily derived for spherically-symmetric winds 
fe.g.. lProgal ll999\ This comparison shows the achievable disk-wind mass loss 
rates differ only from the well-known one-dimensional analytic values by a factor 
(of geometric origin) o f orde r unity. A very similar conclusion was reached by 
iFeldmeier Sz Shlosmanl ( 19991 ) who compared th eir analyt ic LD disk wind models 



the two-dimensional numerical simulations of Ipso dah and found an overall 
good agreement in the streamline shape, tilt angle, and mass-loss rate. These 
agreement between numerical models and analytic ones allow us to generalize 
beyond the limite d set of mod els for which numerical results already exist. 

For example, iDrew Sz Proga (2000) compared mass loss rates predicted by 
the models with observational constraints (Figure 2). They concluded that either 
mass accretion rates in high-state CVs are higher than presently thought by a 
factor of 2-3 or that radiation pressure alone is not quite sufficient to drive the 
observed hypersonic flows. The difficulty in accounting for the mass loss rate in 
a pure LD disk wind model for CVs is simply a reflection of the fact that the CV 
luminosities just barely satisfy the basic requir ement, i.e., Luv ^ 7 x 10~^LEdd- 
A very similar conclusion has been reached bv lMauche &: Ravmondl llOOO) who 
analyzed observations of OY Carinae taken by the Extreme Ultraviolet Explorer. 
iMauche Ravmondl argued that line driving alone falls an order of magnitude 
short of driving the observed mass-loss rate. 

Generally, one can argue that in all accretion disks, with the UV lumi nosity 
L\jy ^ a few 10~^LEdd mass outflows have been observed (jProgall200^ ). For 
example, accretion disks around: massive black holes, WDs (as in AGN and CVs 
with Luv ^ O.OOlLEdd) and low mass YSOs (as in FU Ori stars with Luv ^ 
a few X O.OlLsdd) show powerful fast winds. Systems that have too low UV 
luminosities to drive a wind include accretion disks around neutron stars and 
low mass black holes as in low mass X-ray binaries and galactic black holes. 
These systems indeed do not show outflows similar to those observed in CVs, 
AGN and FU Ori. 

Comparing the wind properties (e.g., and geometry) inferred from ob- 
servations with the wind properties predicted by models for given systems pa- 
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Figure 2. Model mass loss rates for ra diation-driven d i sk win ds as a func- 
tion of r where T = L/LEdd [Figure 1 in lDrew & Progal l|200flD ]. The solid- 
line curve is the 1-dimensional mass loss rate, derived analytically, in the case 
that Mmax — 4000 and a = 0.4. Also shown are numerical disk-wind mass 
loss rates f or the same Mm for a = 0.4 (crosses), 0.6 (circles), 0.8 (aster- 
isks). (See l)Progalll999j) for further details). The vertical Hue superimposed 
is drawn at F corresponding to a mass accretion rate of 10~® Mq yr~^ onto a 
CO white dwarf of mass 1 Mq - presently mass accretion rates are believed to 
be less than this. The solid horizontal line typifies current M ^(73+ estimates, 
while the dashed line is the lower limiting mass loss rate construed from wind 
ionization models. 



rameters is one way of testing the models. It would be very instructive to com- 
pute synthetic spectra based on the models and compare these with observed 
spectra. Such synthetic hne profiles were computed by Proga ct aL ( 2002. ) using 
a generalized version of the Sobolev approximation. In this study, the attention 
was restricted to the case of a representative UV transition of a light ion such 
as C IV or Si IV. The assumed abundance and atomic data were appropriate to 
the C IVA1549 treated as sin glet. 

Generally, iProgal (l2002''l found that the two main wind components (slow 
wind and fast stream) produce distinct spectral features. The fast stream pro- 
duces profiles which show features consistent with observations. These include 
the appearance of the classical P-Cygni shape for a range of inclinations, the 
location of the maximum depth of the absorption component at velocities less 
than the terminal velocity, and the transition from net absorption to net emission 
with increasing inclination. However the model profiles have too little absorption 
or emission equivalent width compared to observed profiles. This quantitative 
difference between the models and observations is not a surprise because, as we 
discussed above, the LD wind models predict a mass loss rate, mostly due to 
the fast stream, that is lower than the rate required by the observations. 

A key parameter shaping the total hne profile - made up of both scattered 
emission and absorption - is the ratio of the expansion velocity to the rotational 
velocity. This ratio is a cause of some differences in the line profile between 
models with or without radiation from the central star. For models with the WD 
radiation switched on (the bottom panels on Figure 3), the winds are less bipolar 
and so the rotational velocity decreases along the wind streamlines faster than 
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Figure 3. Maps of poloidal velocity (the left hand side column of panels) and 
line profiles (2-6 columns of panels) for two representative models of LD disk 
winds. The line profiles are show for a range of of inclination angle, i (see top 
right corner of each panel for the value of j). The top panels are for the wind 
model with x = and Ma = tt x 10"^ M© yr'^ [B2 model in lProsral l|2003bD ]. 
The bottom panels show results f or the m odel with Ma = tt x 10~^ Mq yi:~^ 
and a; = 1 [C2 model in iProgal l|2nn3bll ]. The zero velocity corresponding 
to the line center is indicated by the vertical line in the panels shown line 
profiles. Note the difference in the velocity and flux ranges in the planes for 
i = 85° (the right hand side column). 



for models with the WD radiation switched off (the top panels in Figure 3). This 
simple change in the wind geometry reduces the rotational velocity of the flow 
compared to the expansion velocity. The relatively higher expansion velocity 
has an important consequence on the scattered emission: for x ^ 1 at high 
inclination (i ^ 60°), the red component of the scattered emission becomes 
stronger - stronger than the blue component of the emission and stronger than 
the blueshifted absorption so it is strong enough for the total line to have a 
P Cygni profile (e.g., Figure 3i) 

In the follow up paper, we found that the inclusion of a disk wind at larger 
radii changes qualitati vely arid quan titatively the line profiles predicted by the 
LD disk wind model (iProgal 1200313) ■ The models computed on a small grid - 
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such as those in lPSD 99l . where the outer radius equals 10 WD radii - suffice to 
calculate the gross properties of the disk wind. For example, the radial range of 
10 WD radii suffice to calculate and the fast part of the wind, which are both 
associated with the outflow from the innermost disk. However, such calculations 
do not capture the entire region where lines are formed. As a result, they 
underpredict the line absorption and to a lesser extent the scattered emission. 
The simulations on the small grid predicted a double- hum ped structure near the 
line center for intermediate inclinations fe.g.. lProgall200^ . This structure is due 
to a non-negligible red-shifted absorption that is formed in the slow wind where 
the rot ational velocity dominates over expansion velocity. 

In iProgal ( 2003bl l. we showed that by taking into account the downstream 
part of the same slow wind one is able to increase significantly the central ab- 
sorption. As a result, the double-humped structure is reshaped to a more typical 
broad trough. We emphasize that all improvements in the shape as well as the 
strength of the absorption were achieved without changing the gross properties 
of the wind. In particular, our new models do not predict a higher mass-loss 
rate than the previous models. The changes in the line profiles are mainly 
caused by the fact that the ratio between the rotational and poloidal velocity 
decreases downstream. Overall, one finds that the wind-formed line profiles 
seen at ultraviolet wavelengths cannot originate in a flow where rotation and 
poloidal expansion are comparable. The UV lines must trace gas that expands 
substantially faster than it rotates. 

The main discrepancy between the predicted line profiles and the observed 
ones is in the line emission. Specifically, the model cannot produce the red- 
shifted emission as strong as that seen, for example, in the Civ profile of many 
systems with intermediate inclinations (see below). However, this shortcom- 
ing is not a great surprise - thi s has been a problem for a while (see e.g., 
iMauche. Lee. Sz Kallman 1997; Ko et al.lll996l . and discussion below). 

A systematic comparison between predicted line profiles and observations 
for many systems is a crucial test the idea that the LD disk wind model can 
work for CV winds in the sense that it can reproduce the observed line pro- 
files for model parameters (e.g., L) suitable to CVs. Preliminary results from 
limited survey of dynamical models and their predictions are promising. For 
example. Figure 4 presents a comparison between profiles de rived f r om the 
model with Ma = vr x 10"^ Mgyr^^ x = 0.25 [model E2 in iProgal (|20033) 
for which M„ = 8 x 10~^ M(T^yr~ ^] and observations of the Civ 1549A transi- 
tion of IX Vel ( Hartlev et al.ll200^ . [To show how much line emission is required, 
only the absorption cornponen t is plotted.] The observed i for this system is 60° 
( Beuermann &: Thomaslfl990t ). 

Figure 4 shows that the model proffies well reproduce the blue-shifted ab- 
sorption despite a relatively low M^. Thus the gap between the kinematics of 
the LD wind models and reality is narrower than comparison between observed 
and theoretical M^'s would indicate. The gap is only narrowed but not bridged 
yet because the mass fluxes required to match the o bserved spectra at least of 

IX Vel are somewhat higher than those observed. As lDrew Sz Progal (|2P00) dis- 
cussed, the luminosity of the system requires a mass accretion rate of at most 
10~* M0yr~^, whereas the line profiles require Ma = vr X 10-s M© yr-i and 

X = 0.25, yielding a system luminosity higher than the observed one by a factor 
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Figure 4. A map of poloidal velocity and line profiles. The right hand side 
panel, compares profiles derived from a wind model for i — 50° and 60° (thick 
solid and thick dashed line, respectively) and observations of t he CiV 1549A 
in th e spectrum of the brightest nova-like variable, IX Vel llHartlev et al.l 
|2002|) . The synthesized lines show the absorption component without the 
contribution from the scattered emission (see the main text). 



of ~ 4. However, the main point here is that this discrepancy is much smaller 
than it used to be (i.e., it was more than 1 order of magnitude) and there is a 
good chance that it can be r educed still further. 

For example, in IProgal (2QQ23) we have computed many models, changing 
various model parameters suc h as Mp and th e parameters of the force multiplier, 
a and Mmax- In general, in IProgal ( 20n3lj ) found that there is a degeneracy 
in the model parameters as far as line profiles are concerned fPSD 98 found 
an analogous degeneracy for the wind properties). For example, models with 
slightly different parameters - such as Ma and a - produce similar line profiles for 
different i. Additionally, a model with with Ma = 10~^ Mq yr"-*^ and a = 0.674 
predicts very similar line profiles to model with Ma = vr x 10^^ Mq yr^^ and 
a = 0.6 for the same i. Finally, we find that the product (1 -|-a;)-LDMjnax) not its 
individual factors, appears to be a fundamental parameter determining the line 
profiles (i.e., their width and depth) for the parameter range applicable to CVs. 
This has an important implication for LD models: to obtain a theoretical fit as 
good as shown in Figure 4 for a fixed i, one needs (l-|-a;)LDMmax ~ 1-3 x 10^ Lq 
rather than specifically x = 0.25, Mmax = 4400, and Ld = 23.4 Lq as for the 
model shown in Figure 4. 

The above mentioned successes of LD disk wind models are encouraging but 
there are also problems. For example, if radiation pressure powers the mass loss 
from CVs, t he wind mass-loss ra te should increase with increasing system lumi- 
nosity (e.g. JPSD 98l : lProgJl999l ).). However recent HST observations of IX Vel 
and V3885 Sgr showed that the wind spec tral features do not correlate with 
the system luminosity ( Hartlev et al.ll2002t ). If confirmed, these observations 
seriously challenge the pure LD disk wind scenario. 
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2.3. Limitations of LD models and future work 

We conclude this section with a reminder of the hmitations of the present LD 
disk wind models. Many of the details of the models depend on details of 
the assumptions about the disk and microphysics of the wind. For example, 
the assumed form for the background radiation field inevitably plays a role in 
the comparison of observed profiles with synthetic profiles based on any model, 
kinematic or dynamical. However for line fitting based on LD wind models, the 
adopted disk and WD radiation fields are even more important because they 
determine all wind properties except the initial Keplerian component of motion. 
Current models adopt the dependence of the d isk radiation o n radius according 
to the standard steady state disk model (e.g.. lPringlelll98lh . This assumption 
is a good starting point but it should be remembered that actual disks may 
only be very crudely described by simple theory or they may yield unexpected 
features such as chromospheric emission. In fact, spectra l synthesis models of ac- 
cretion disk p hotospheres fe.g.. iLinnell Sz HubenvllToQ^ : IWade fc Hubeiivlll998l : 
IWade Sz Orosa..l999 ) have typically failed to adequately reproduce o bserved en- 
ergy distributions where direct comparisons have been made (see e.g.. lLong et al.l 
Il994 ). There has been a similar lack of success in past calculations of accretion 
disk line emission. Specifically, the accounting for both the strengths of and 
the flux ratios among various emission lines in a large ensemble of observe d 
CVs has not been compelling ( Mauche. Lee. Sz Kallmanlll997l : IKo et al.lll996t ). 



These models either miss a crucial physical component or employ an inappro- 
priate physical assumption which afi'ects the predicted line emission. It is for 
this reason that, in comparing model wind profiles with observation, we concern 
ourselves more with the absorption component than with the emission. 

In line profile calculations, Proga (2002) and Proga (2003b) assumed that 
the ionization fraction of scattering species is constant (^ion is set to unity). A 
proper allowance for a variation of ^ion with position to values less than one 
can only serve to weaken the overall line profile. [Examples of plots of this 
potentially very ma r ked p ositio nal dependence may b e found in the work of 
IShlosman fc Vitellol (Il99.'j 'l and iLong fc Kniggd (|2002l 'l]. Therefore, it is im- 
portant to carry out time-dependent calculations of the wind photionization 
structure. It is also important to perform simulations in fully three dimensions 
to explore nonaxisymmetric effects. The two photoionization and three dimen- 
sional effects are most likely coupled (e.g., by the wind density) and it would be 
essential to study them self-consistently. 



3. MHD Models 

One of the reasons for considering magnetic fields as an explanation for winds 
from accretion disks is the fact that magnetic fields are very likely crucial for 
the existence of all accretion disks. The magnetorotational instability (MRI) has 
been shown to be a very robust and universal mechanism to p roduce turbulence 
and the tra nsport of angular momentum in disks at all radii tealbus Sz HawlevI 
119911 . [ 19981 1. It is therefore likely that magnetic fields control mass accretion 
inside the disk and play a role in producing a disk wind. 

In fact, magnetically driven winds from disks are t he favored explana- 
tion for the outfiows in many astrophysical environments. iBlandford Sz Pavnd 
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( 19821 ) (see also iPelletier fc Pudritz I Il99^ 'l showed that the centrifugal force 



can drive a wind from the disk if the poloidal component of the magnetic field, 
Bp makes an angle of > 30° with respect to the normal to the disk surface. 
Generally, centrifugally-driven MHD disk winds (magnetocentrifugal winds for 
short) require the presence of a sufficiently strong, large-scale, ordered mag- 
netic field threading the disk with a poloidal comp onent at least comparable 
to the toroidal magnetic field, \B(i^/Bp\ ^ 1 (e.g., ICannizz~ fc Pudritd IT989I : 
IPelletier fc Pudritz l ll99'3). Several groups have studied n umerically ax isymmet- 

ric winds using the Blan d ford fc Payne mechanism fe.g.. | Usvugova"et al...l995l: 

Ouved fc Pudritzll997al lbl: lKrasnoDolskv. Li. fc Blandfordll999l : lKato. Kudoh fc Shibatal 
2OO2I ). An important feature of magnetocentrifugal winds is that they require 



some assistance to flow freely and steadily from the surface of the di sk, to pass 
through a slow magnetosonic surface (e.s.. iBlandford Sz Pavnelll982l ). The nu- 
merical studies mentioned above do not resolve the vertical structure of the 
disk but treat it as a boundary surface through which mass is loaded on to the 
magnetic field lines at a specified rate. 

Winds from disks can also driven by the magnetic pressure. In particular, 
the toroidal magnetic field can quickly builds up due to the differential rotation 
of the disk so that \B^/Bp\ » 1. In such a case, the magnetic pressure of the 

toroidal field can give rise to a self-starting wiiid fe.g.. lUc hida &: Shibata 19851: 

Pudritz fc Normanll98(^ : IStone fc Normanll994l : IContoDoulosii,199,^ : .Kudoh fc Shibatal 



1991 ). To produce a steady outflow driven by the magnetic pressure a steady 



supply of advected toroidal magnetic flu x at the wind base is needed, oth erwise 
the outflow is likely to be transient fe.g.. lKonig]||l993l : IContopouloslll995l ). It is 
still not clear whether the differential rotation of the disk can produce such a 
supply of the toroidal magnetic flux to match the escape of magne tic flux in the 
wind and even if it does whether such a system will be stable fe.g.. IContopoulo^ 
Il99 5: Ouyed & Pudritz 1997b, a nd ref e rences therein). 

To our best knowledge, in IProgal ( 2003al ). we were first to report on nu- 



merical simulations of the two-dimensional, time-dependent MHD structure of 
LD winds from luminous accretion disks initially threaded by a purely axial 
magnetic field. We developed self-consistent models of such winds and applied 
them to winds from CVs. Our models require less free parameters than previous 
MHD models . In par ticular, the model predicts the mass loss rate. 

In lProgal ( 2003al ). we used ideal MHD to compute the evolution of Keplerian 



disks, varying the magnetic field strength and Ld, L^, or both. We found that 
the magnetic field very quickly starts deviating from purely axial due to MRI. 
This leads to fast growth of the toroidal magnetic field as field lines wind up due 
to the disk rotation. As a result the toroidal field dominates over the poloidal 
field above the disk and the gradient of the former drives a slow and dense disk 
outfiow, which conserves specific angular momentum of fluid. 

Our LD-MHD simulations also showed that depending on the strength of the 
magnetic field relative to L the disk wind can be LD or MHD-driven. For very 
weak magnetic fields, similarity to the LD wind, the wind consists of a dense, 
slow outflow that is bounded on the polar side by a high-velocity stream. The 
mass-loss rate is mostly due to the fast stream. As the magnetic field strength 
increases first the slow part of the flow is affected, namely it becomes even denser 
and slightly faster and begins to dominate the mass-loss rate. In very strong 
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magnetic field or pure MHD cases, the wind consists of only a dense, slow outflow 
without the presence of the distinctive fast stream so typical of pure LD winds. 
Our simulations indicate that winds launched by the magnetic fields are likely 
to remain dominated by the fields downstream because of their relatively high 
densities. Line driving may not be able to change a dense MHD wind because 
the line force strongly decreases with increasing density. 

The increase of the mass loss rate due to the MHD effects is a welcome 
development in modeling CV winds. However, it is unclear whether LD-MHD 
models can resolve the problem of too low M„. As we discussed in section 2, 
to explain CV winds we need a model that predicts not only a higher AI^ but 
also M„ must be mostly due to a fast wind not a dense slow rotating wind as 
we found in our LD-MHD models. 

3.1. Limitations of LD-MHD models and future work 

The most important limitation of LD-MHD simulations is an inadequate spatial 
resolution for modeling the MRI inside the disk. These are only preliminary 
simulations we aimed to examine the parameter space of the models that will 
define the major trends in disk wind behavior. Therefore the priority has been so 
far to set up the simulations in such a way that the base of the wind is relatively 
stable and corresponds to a steady state accretion disk. 

The fact that LD-MHD results strongly depend on the magnetic field points 
to a need to explore different configurations for the initial magnetic field and to 
move from two-dimensional axisymmetric simulations to fully three-dimensional 
simulations. It is important to follow the long-time evolution of the flow. There- 
fore three-dimensional simulations are required as there exist no self-sustained 
axisymmetric dynamos. Thus, contrary to the stellar winds, simulations of winds 
from magnetized disks - with or without radiation pressure - should include the 
disks themselves, not just the disk photosphere, and should be performed in 
three dimensions. 



4. Conclusion 

Our main conclusion is that, despite some problems, line-driving alone is still the 
most plausible mechanism for driving the CV winds. Preliminary results from 
LD-MHD wind models confirm that magnetic driving is likely an important 
element of the wind dynamics. However, magnetic driving does not seem to 
be necessary to produce a wind. The most important issues which need to be 
addressed by future dynamical models, regardless of driving mechanism, are the 
effects of the position-dependent photoionization and the dynamical effects in 
three dimensions. 
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